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Abstract 

We analyze the statistics of electrostatic energies (and their differ- 
ences) for a quantum dot system composed of a finite number K of 
electron islands (metallic grains) with random capacitance-inductance 
matrix C, for which the total charge is discrete, Q = Ne (where 
e is the charge of an electron and N is an integer). The analy- 
sis is based on a generalized charging model, where the electrons 
are distributed among the grains such that the electrostatic energy 
E(N) is minimal. Its second difference (inverse compressibility) xn = 
E{N + 1) — 2E(N) + E(N — 1) represents the spacing between ad- 
jacent Coulomb blockade peaks appearing when the conductance of 
the quantum dot is plotted against gate voltage. The statistics of 
this quantity has been the focus of experimental and theoretical in- 
vestigations during the last two decades. We provide an algorithm for 
calculating the distribution function corresponding to xn and show 
that this function is piecewise polynomial. 

1 Introduction 

The physics exposed in the addition spectra of quantum dots is rather rich, 
and hence its investigation is at the focus of both experimental and theoreti- 
cal studies. After the origin of Coulomb-blockade peaks has been elucidated, 
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investigation is directed toward more subtle questions such as their heights, 
widths, and spacings. The underlying physics is related to the ground-state 
energy, chemical potential, and inverse compressibility of quantum dots com- 
posed of a few metallic electron islands coupled capacitively and inductively 
to each other. 

The present work concentrates on the distribution of spacings between 
Coulomb-blockade peaks in large semiconductor quantum dots. In particu- 
lar, we are interested in fluctuations of these quantities with the number N 
of electrons on the dot. The main problem can be stated as follows: Accord- 
ing to the simplest picture (charging model), in which the quantum dot is 
regarded as a single-electron island whose coupling with the leads is through 
its capacitance C, the total potential energy of a quantum dot with iV elec- 
trons and charge Q = Ne is Q 2 /2C — V g Q, where V g is the gate voltage and 
e is the electron charge. The position of the iV-th Coulomb-blockade peak 
occurs at a gate voltage V g = ^f-. This peak position is then a linear func- 
tion of N, and therefore the spacing between two adjacent peaks should be a 
constant e 2 /C , independent of N. This is not always confirmed experimen- 
tally. The situation is even more intriguing if the quantum dot is large and 
might contain more than a single electron puddle. As indicated in a series 
of recent experiments [T], the spacing between adjacent Coulomb blockade 
peaks occasionally vanishes; namely, Coulomb blockade peaks tend to bunch. 
The problem is therefore to explain why the results predicted from a simple 
charging model deviate substantially from the experimental observation. 

In j2], a generalized charging model has been tested, where it is assumed 
that the large dot used in the experiments jlj could be divided into a set of 
potential wells (metallic grains) with random capacitances and random mu- 
tual inductances. This casts the question of Coulomb blockade peak spacing 
distribution into the problem of elucidating the statistics of the addition 
spectrum of a relatively simple physical system. It consists of K metallic 
grains (or capacitors), such that the number of electrons on the z-th grain 
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is rij (z = 1, 2, . . . , K), the total number of electrons being N. The charging 
model for such a system (at zero temperature) is based on the assumption 
that the distribution of electrons among the grains is determined by requir- 
ing that the electrostatic energy E(N) of a dot containing N electrons is 
minimal. It is useful at this point to recall the basic facts pertaining to the 
energy of the electrostatic field of conductors [3]. The electrostatic energy 
of the system is a bilinear form in the numbers n*. This form is given by a 
K x K matrix W = \C~ 1 . Here C (matrix elements c^-; i,j = 1,2. ..if) 
is a positive-definite symmetric matrix of capacitance and inductance coef- 
ficients. Physically, the matrix C has positive diagonal entries and negative 
(more precisely, non-positive) non-diagonal entries [3], 

Cij = Cji, cu > 0, Cij < (i^ j). (1) 

On the other hand, all the elements of C^ 1 are non-negative. More precisely, 
with Wij] i,j = 1, 2...K the elements of the matrix W = |C _1 , one has, 

wu > 0, w tj > (i^j). (2) 

The off-diagonal entries c^-, i ^ j, decay as an inverse power of the distance 
between the grains, while the diagonal entries cu are proportional to the 
geometrical size of the grains. The notion of randomness enters when we 
recall that experimentally, the sizes of the grains, as well as the distances 
between them, are random quantities. This means that the elements of the 
matrix C are random numbers (subject, of course, to the required symmetries 
((H)). The spacing between Coulomb blockade peaks is equal to the second 
difference of the ground state energy. In other words, the distribution of 
spacing peaks is determined by the statistics of the inverse compressibility, 

XN = E(N + 1)-2E(N) + E(N-1). (3) 

When two Coulomb blockade peaks coincide, the second difference xn van- 
ishes. Note that, on the average (and on a large scale), the energy E(N) grows 
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quadratically with N. Therefore, one would expect the second difference to 
stay finite and N independent. However, there is no simple relation such as 
E(N) = aN + bN 2 . The deviation of E(N) from exact quadratic behavior 
makes its second difference xn non-constant, and a fluctuating quantity. It 
is precisely these fluctuations which we intend to study. As we shall see, the 
fact that electron charge is quantized makes this task non-trivial. 

Having explained the physical motivation, we then pose the mathematical 
problem: what is the distribution p(x) of inverse compressibility for a given 
system of metallic grains with random capacitive matrix C? As a crude 
approximation it was assumed in [3] that the metallic grains are indeed very 
far apart, and the matrix C is nearly diagonal, its K diagonal elements 
(capacitances) being random numbers. The energy of the system in the 
diagonal case is given by 

K 1 

E(N) = min ^ — n?, (subject to J2? =1 rii = N). (4) 
i=i ICl 

The minimum in (J3J) is taken over all possible partitions (rii)f =l of N. It 
was first proved that the minimum problem (j3J) has the following convenient 
feature: If n 1; n 2 , . . . , Uk are the argument values bringing E(N) to its min- 
imal value for some N, then the minimum for N + 1 is obtained by retaining 
all rij's, except for one which is increased by 1. This allowed an exact deter- 
mination of the distribution function according to which the sequence xn is 
distributed. For a random set of capacitors (ci, c 2 , ck - random numbers 
with probability distribution P(ci, c%, ck)), the distribution of the inverse 
compressibility F(x) was calculated in 0. 

Our next goal is to study this problem for general positive-definite ma- 
trices C. The problem turns out to be quite harder. To begin with, it is 
no longer true that the optimal solution for N + 1 is obtained in a simple 
manner from that for N. That is, for each N we need to re-distribute the N 
electrons between the grains, and it may well happen that, although the n^s 
grow in general with N, some of them will actually decrease infinitely of- 
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ten as N increases by 1 each time. Namely, there will exist infinitely many 
values of iV for which the optimal value of some n» decreases as N grows 
to JV + 1. Thus, the problem entails new behavioral patterns with respect to 
the diagonal case. 

Our main result in this paper is an algorithm for calculating the distribu- 
tion function corresponding to (xn)- Moreover, we show that this function 
is piecewise polynomial. We state the result in Section 2. Section 3 is a short 
digression, discussing a few notions which arise in the proof. The proof of 
the main theorem is given in Section 4. 



2 The Main Results 



Mathematically, our problem is as follows. Let C = (cij)fj =1 be a positive- 
definite matrix, with positive diagonal elements and non-positive off-diagonal 
elements. Assume that the sum of elements in every row of C is positive and 
that all entries of \C~ l = W = {wij)f^ =1 are non-negative. Put 

{k k \ 

^2 WijUiUj : m G Z+, J^n; = N > , iV G N, (5) 
i,j=i i=i J 

where N is the set of positive integers and Z + = N U {0}. We want to 
understand the statistical behavior of the sequence E(N), and in particular 
that of the second difference sequence 

X n = E(N + 1)-2E(N) + E(N-1). (6) 

To formulate our main result, we need a few definitions and notations. 

Definition 2.1 Let (x n )~ =1 be a sequence of real numbers and F a distri- 
bution function. The sequence (x n ) is asymptotically F- distributed if 

|{l<n<M:x w <x}| _^ 

M M^co 
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for every continuity point x of F (where \S\ denotes the cardinality of a finite 
set S). 



The definition almost coincides with [5J p. 53, Def. 7.1], except that there 
the sequence (x n ) is considered only modulo 1. Note that a sequence need 
not be asymptotically F-distributed for some F, as the following example 
shows. 

Example 2.1 The sequence of numbers 



is not asymptotically F -distributed for any F. 

A stronger notion is obtained when we require not only long initial block 
of the sequence to behave approximately according to F, but rather require 
any long block to behave so. This leads to the following definition (5, p. 40, 
Def. 5.1] and p.200, Def. 3.2].) 

Definition 2.2 In the setup of Definition 12.11 (x n ) is asymptotically well 
F-distributed if 



for every continuity point x of F. 

The following example demonstrates that the property of asymptotic 
well F-distribution is indeed strictly stronger than that of asymptotic F- 
distribution. 

Example 2.2 The sequence of numbers 



0, 0,...,0, 1,... 



1,... 




\{L < n < M : x n < x}\ 
M-L 



M— L— too 



F{x) 







1 , 0,0, 1,1,0,0,0,1,1,1,... 




1 



1 



2 



2 



3 



3 
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is asymptotically F -distributed, where F is the distribution function 

ro, x<o, 

F(x)=l I, 0<x<l, 
[l, x>l, 

but it is not asymptotically well F -distributed. 

Definition 2.3 A function g : R — > R is piecewise polynomial if there exist 
intervals (finite or infinite) Ij C R and polynomials Qj, 1 < j < m, such 
that 

g(x) = Qj(x), x e Ij, 1 < j < m. 
The degree of g is max cleg Qj. 

l<j'<m 

Returning to our problem, let 6j = 52j=i c «i' ^ — * — the row sums 

of the matrix C. Since the matrix C is random, in the generic case the 
numbers b±, 62, •••) bx are linearly independent over the rationals. (That is, 
considered as vectors in the space R over the field of rational numbers Q, 
they are independent.) 

Now we can formulate our main result. 

Theorem 2.1 Let C be a positive-definite symmetric matrix, with positive 
row sums 61, b 2 , ...,bK and let W = \C~ l . Suppose that bi, 6 2 , bx are lin- 
early independent over the rationals. Then the sequence (xn)n=i °f ^ ne sec " 
ond differences, defined via © and (0), is asymptotically well F-distributed, 
where F is a continuous piecewise polynomial function of degree at most 
K — 1, which can be effectively computed. 

As mentioned in the introduction, a phenomenon which occurs in the gen- 
eral case dealt with here, but not in the special case of diagonal matrices C, 
is that, as we pass from N to N + 1, there may be re-distribution of the n^s 
in the optimal solution. The following example is to that effect. 
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E(N) = — min (En-, + 2riin 2 + Ariin s + 5n 9 + An^ri?. + Arioi. 

16 n 1 +n 2 +n 3 =N,n i >0 

Using the techniques in the beginning of Section 4, it is easy to verify that 
the optimal values of ni,n2,n 3 are given by 



A N N ^ AT = 0(mod3), 
(n u n 2 ,n 3 ) = { (^1,^1, M?), iV^l(mod3), 



iV = 2(mod3). 

Thus, for any non-negative integer K , when passing from N = 3k + 1 to 
N = 3k + 2, the value of n 3 at the optimal point decreases from k + 1 to k. 

3 Uniform Distribution Modulo 1 

In this section we briefly discuss the notion of uniform distribution modulo 1 
and recall a few related results, which will be needed in the proof of the 
Theorem 2.1. 

Definition 3.1 A sequence (x n )^_ 1 of real numbers is uniformly distributed 
modulo 1 if 

\{l<„<N:a< {x „}<b}\ _ <„<6<1, 

where {t} is the fractional part of a real number t (0 p.l, Def. 1.1]). 

In terms of Definition 12.11 (x n ) is uniformly distributed modulo 1 if and 
only if the sequence ({x n }) of fractional parts is F-distributed, where F is 
the distribution function of the uniform distribution on [0, 1]: 

x < 0, 

F(x) = { x, < x < 1, (7) 
x > 1. 



8 



The notion of uniform distribution modulo 1 has a multi-dimensional 
analogue. A sequence (x n )^ ( L 1 in R s is uniformly distributed modulo 1 in R s 
if 

|{1 < n < N : a < |x n | < b}| -A-,, 

11 = = ~ 1 1 — > T[(h - a *)' < a < b < 1, 

iV TV— »oo -*- A 

i=l 

where = (0, 0, . . . , 0) G R s , a = (ai, a 2 , . . . , a s ), and so forth, and inequal- 
ities between vectors in R s are to be understood component-wise (0 p. 47, 
Def. 6.1]). 

The notion of uniform distribution modulo 1, both in the 1-dimensional 
and the multi-dimensional cases, has a stronger version, whereby the required 
property holds not only along initial blocks of the sequence, but along any 
blocks of larger and larger lengths (0 p. 40, Def. 5.1]. A sequence satisfying 
this stronger property is well distributed modulo 1. Obviously, well distribu- 
tion modulo 1 is equivalent in the 1-dimensional case to F-distribution of the 
sequence of fractional parts for the function F given by (j7J) . A basic example 
of a sequence which is uniformly distributed modulo 1 is (na)™ =1 , where a 
is an arbitrary irrational p.8, Def 2.1]. In the mult i- dimensional case, the 
sequence (nai,na 2 , ■ ■ ■ , , n.a s )^i 1 is uniformly distributed modulo 1 in R s if 
and only if the numbers 1, a\, ■ ■ ■ , at s are linearly independent over Q [5, 
pp. 48-49]. Moreover, for these sequence, well distribution is equivalent to 
uniform distribution. 

Recall that the density of a set A C N is given by 

D(A) = km J — — - 

V ' M->oo M 

if the limits exists. If, moreover, the limit 



BD(A) = lim 



\An{L,M}\ 



M-L-*oc M 

exists, then it is called the Banach density of A. 

We can rephrase the definition of uniform distribution modulo 1 using the 
notion of density of a set. Namely, (x n )^ ( L 1 is uniformly distributed modulo 1 
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if for every interval I C [0, 1) we have 



D({n : {x n } e I}) = \I 



(8) 



where |/| denotes the length of /. Similarly, (x n )^ =1 is well distributed mod- 
ulo 1 if (jHJ) continues to hold when the density of the left-hand side is replaced 
by Banach density. 

4 Proof of Theorem 2.1 

To avoid complicated notation, shall prove in Theorem 2.1 only that (xjv)^ =1 
is asymptotically F-distributed, and not that it is asymptotically well F- 
distributed. As will be seen in the proof, our result depends on the fact that 
the sequence ({Nb\}, {Nb 2 }, {Nba-i}) is uniformly distributed modulo 1 
in H K ~ 1 . Since this sequence is actually well distributed modulo 1, the same 
proof shows that (xn)n=i * s actually well F-distributed. 
Along with the sequence E(N) from it is very useful to consider the 
sequence E%(N), defined by 



Obviously, Ei(N) < E(N) for each N. We shall refer to the minimum 
problems on the right hand side of (0) and of © as the constrained problem 
and the unconstrained problem, respectively. 

Denote by e the column K- vector with all entries 1. 

Lemma 4.1 The unique minimum of the unconstrained problem is 



Ei(N) = mini 



e R, J2 x i = N \i N eN. (9) 



i=l 



N 



■Ce. 




and the corresponding unconstrained minimum is Ei(N) 
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Due to our assumption regarding the positivity of the row sums of C, all 
components of xo are positive. Multiplying all entries of C by any constant 
c > we obtain an equivalent problem. Taking c = (5^ij=i c ij)~\ we 
shall henceforth assume that X]fj=i c ij = 1- m particular, denoting 
b = (&i, 6 2) ■■■bK) t , we have 

x = iVCe = Nh (10) 



and 



AT 2 

^iW = — • (11) 



Proof of Lemma 14.11 Let x ^ xo be any feasible solution of the uncon- 

N 

strained problem. Putting a = — ^ and y = x — x , we obtain 

x'lfx = (x + yYW(x + y) = x&Wx + 2x£Wy + y*Wy 
= x*^x + al*CC- x y + y'Wy 

= x* VTx + ol*y + y l Wy = x^Wxo + y'Wy > x^Wxq . 

Consequently, 

N N N 2 
E 1 (N) = —f e'CW— fP- Ce 



Lemma 4.2 E 1 (N) > E(N) - J2f,j=i Wi i f or 

every N. 

Proof Let x = (xoi,x 02 , . . . ,x k) be the minimum point of the uncon- 
strained problem. Let r = Y^i=ii x oi} be the sum of fractional parts of all 
coordinates of x . Obviously, r is an integer, < r < K. Let i 1: i 2 , . . . , %k be 
all integers between 1 and K, ordered so that {xo.ii } < {^0,12 }<•••< { x a,i K } 
(where ties are resolved arbitrarily). Consider the vector n = (m, n 2 , . . . , uk) 
defined by 

[X 0t i], i = h,i 2 , ■ ■ ■ ,lK-r, 

[xo t i] + 1, otherwise . 
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As mentioned in Lemma 14. H all Xq/s are positive, and hence n is a feasible 
solution of the constrained problem. Set y = n — xo. Since all coordinates 
of y lie in the interval (—1, 1), as in the proof of Lemma f4. II we have 

K 

E(N) < n'Wn = x^Wtcq + y l Wy < E X {N) + ^ Wij (12) 
which proves the lemma. 

Lemma 4.3 There exists an effective constant A = A(C) such that, for 
every N , the distance between the solution of the constrained problem and 
that of the unconstrained problem does not exceed A. 

Proof Write W = P~ 1 DP, where P is orthogonal and D diagonal. Let 
M be an upper bound on the eigenvalues of C (for example, the L°°-norm 
Ef=iM of C). Then M" 1 is a lower bound for the eigenvalues 
of W, namely for the diagonal entries of D. Let F be the diagonal matrix 
with positive diagonal entries and F 2 = D. Obviously, H-Fz^ > Af _1 / 2 ||z|| 2 
for every z G R . Then for every y G R^ we have 

y t Wy = y t p- l FFPy=\\FPy\\l>\\Py\\l/M=\\y\\l/M. 

Now let x and n = x + y be minimum points of the unconstrained problem 
and of the constrained problem, respectively. Then 

E(N) = n'Wn = x* Wx + y'Wy, 

which implies by Lemma f4. 21 that y f Wy < Y2fj=i w ij- Thus 

K 

l|y||!/M< J2 w v> 

which yields the conclusion of the lemma with 

K 
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Proof of Theorem 2.1. Lemmas I4.1M.3I provide a simple algorithm for 
calculating E(N) for each iV in constant time. Namely, we find the point 
xo yielding the optimal E±(N) according to Lemma f4.H calculate the value 
of i^Wn for all integral points n, with coordinate sum N, within distance 
A from xo, and take the best of them. If the optimal point turns out to 
be n = xo + y, we shall refer to y as the correction vector. We have y = 
n — xo = I — {xo}, where {xo} denotes the vector of fractional parts of 
the coordinates of x and I belongs to some finite effective set L of integer 
vectors. Since the sum of coordinates of the correction vector is always 0, 
the sum of coordinates of I must equal that of {x }. Thus, L consists of all 
integer vectors I, for which the vector I — {x } is of norm not exceeding the 
bound in Lemma 14.31 and its coordinates sum vanishes. To emphasize the 
dependence of L on xo, we shall sometimes write L(xq) instead of L. 

Now when choosing the optimal I out of L, we first notice that, among 
any two candidates l\ and l 2 , the former will be better (or equal) than the 
latter if and only if 

(xq + h- {x })*W(x + h - {x }) < (x + l 2 - {x })'W(x + l 2 - {x }) . 

This inequality is easily seen to be equivalent to 

2(l 2 - IxYWixo} < l\Wl 2 - l\Wh . 

Consequently, / is the optimal choice if and only if 

2(1' - ifWixo} < {IJWV - l*Wl, I'eL. (13) 

To study the second differences 

Xn = E(N + 1) - 2E(N) + E(N - 1), 

we shall write each term on the right-hand side in the form Ei[N + j) + dj 
for an appropriate dj. In fact, as in (fT2|) . denoting by yi,y2 and y 3 the 
correction vectors for N — 1, N and N + 1, respectively, we have: 

X n = E 1 (N + 1)- 2Et (N) + E 1 (N-1) + y* Wy 3 - 2y* Wy 2 + y[W yi . (14) 
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By (HU): 

X n = 1 + ylW2y 3 - 2y\Wy 2 + y\Wy x . (15) 

Let be x , Xq, Xg the points yielding the optimal values of E\(N — 1), Ei(N), 
Ei(N + 1), respectively. In view of (JTUJ): 

x = x + b, x'J = x + 2b. (16) 

For appropriate integer vectors p G £(x ), p' G £(x' ), p" G £(x '): 

yi = P - {x }, y 2 = P ' - {x' }, y 3 = p" - K>- (17) 

The vectors p, p', p" are determined by the system of inequalities: 

2(1 - p) t W^{x } < l l Wl - p'Wp, Jel(xo), 
2(1 - p'YW{x' } < tWl - (p'YWp', I G L(x' ), (18) 

2(Z - p")W{x^} < Z*iy/ - (p"YWp", I G L(xg). 

Due to ()16|) . it is natural to try to rewrite ()18p in terms of {x } without 
referring to {x[,} and {xq}. Divide the K- dimensional torus T , which we 
identify with [0, 1) K , according to the vector b, as follows. 

The i-th coordinate {x' 0i } of {x' Q } may be either {xoi} + h or {xch} + h — 
1, depending on whether {x' 0i } + hi is smaller than 1 or not, respectively. 
Similarly, {x^} may assume one of the three values {x 0i } + 2bi — c, where 
c = 0, 1, 2. Divide the circle T into three disjoint intervals (actually arcs), on 
each of which both {x' 0i } and {x^} assume the same form in terms of {xoi}. 
We have to distinguish between two cases: 
1) If bi < |, write 

T = [0, 1 - 2b t ) U [1 - 2b h 1 - bi) U [1 - k, 1). (19) 

If {xoi} belongs to the first interval on the right-hand, then 

Woi) = i X 0i} + K Wu} = {xoi} + 2&i, 
if it belongs to the second 

{z'oJ = i x oi} + bi, {x'oi} = {x 0i } + 2bi - 1, 
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and if it belongs to the third 

{x' oi } = {x 0i } + bi-l, {xqJ = {x 0i } + 2k - 1. 

2) If bi > |, write 

T = [0, 1 - h) U [1 - k, 2 - 26 ?; ) U [2 - 26i, 1). (20) 

This time, depending on the interval on the right-hand side containing {xoi}, 
we have either 

{x' J = {x 0i } + h, {xq'J = {x j} + 2bi - 1, 



or 



or 



{a4} = {x j} +6.-1, {zqJ = {^oi} + 26i - 1, 



{x' oi } = {x 0i } + bi-l, {xq'J = {x j} + 2k - 2. 

Let In,I i2 ,Ii3 be the intervals on the right-hand side of (fHJj) or (J2UJ), 
depending on whether 6j < | or not, respectively. Denote: 

G{1,2,3}. 

decompose the (K — l)-dimensional torus into a union 
of 3 K ~ 1 disjoint boxes: 



33 3 

o 

V1V2---VK-1- 



t k ~ 1 = |J U - U " 



The information provided by the vector {x } is partly redundant as the 
fact that ^2f = i{xoi} is an integer determines each component in terms of the 
others. To avoid this inconvenience, we shall eliminate, say, {xoa'}- Divide 
T^ 1 into K parts as follows: 



K-x 



Q. s = {{t x ,...,t K -x) eT w :s-l< J2 f ^ s }' s = 0,...,K-l. (21) 



i=i 
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(Thus, tt° = {0}, while all other fi*'s have non-empty interior.) Suppose 
that ({xoi}, {xo,k-i}) € Then: 

{x k} = s - {x 01 } - ... - {x ,a'-i}- (22) 

We need a further subdivision to ensure that, in each cell, both {x' 0K } 
and {x'q K } assume the same form in terms of {xoi}, {£o,i<r-i}- To this end, 
we first split T into three subintervals, similarly to (|19|) and (J2(Jj) . depending 
on whether < ~ or bjc > |, namely 

T = [0, 1 - 2^) U [1 - 2^, 1 - b K ) U [1 - V, 1) (23) 

or 

T = [0, 1 - V) U [1 - b K , 2 - 26*) U [2 - 2b K , 1). (24) 
Let im., Ik2, Iks be the intervals in the splitting. Put: 

K-l 

fi» = {fa, G IT : s - ^ U e W« »7 = 1, 2, 3. (25) 

i=l 

Suppose 6x < §• If the point ({x i}, {xo,k-i}) belongs to flf, then 

Wok) = s ~ Efi^w} + bK, {x'o K } = s- E^NO + 2b K , 
if it belongs to flf, then 

{^dx} = s - YlfJiWi} + {4V} = s - YlfJiWi} + 2b K ~ 1, 

and if it belongs to fig, then 

{x'oiJ = s - E l =T 1 {^oi} + b K - 1, {a&J = s - E^Ni} + 26x - 1. 

If 6# > |, then we similarly find linear expressions for {xq^} an d {^o'ic} i* 1 
terms of the {xch}'s on each fH. 
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Denote: 



Km-m = ^m-m-, n n' VK , < s < K - 1, 1 < Vi < 3. 

Then 

t" 1 U J U 



JsT-1 3 3 

m-VK 

=0 771=1 »?if=l 



forms a decomposition of T into a disjoint union of if • 3 disjoint poly- 
topes. The important property of this decomposition is that, if ({x i}, 
{xo,k-i}) belongs to any cell W mm ... m , the 2K + 1 numbers {x 0K }, Woi}i 
{x' QK }, {x'oi}, • Wok} depend linearly on the first K—l coordinates {xoi}, 
{xq } k-i}- That is 



KJ = {x 0l } + b i -a , i , 1 < % < K - 1, 

KJ = NO + 26, - a?, 1 < % < K - 1, 

K-l 

Wok} = s - ^2i x oi} + b K ~ at' K , 



(26) 



i=l 
K-l 



(27) 



Wok} =8 - ^Wi} + 2J} k 
i=i 

where a[ G {0,1}, a" G {0,1,2}, i = 1,2,..., if. Altogether, there exists 
a linear transformation T : H K ~ 1 — > R^, and for each cell fi* „ there 
exist vectors v, v',v" G R^, such that, denoting x = ({a;oi}> {xo,k-i}), 
we have 

{x } = Tx + v, {x^Tx + v', K} = Tx + v". (28) 

On each cell we may now rewrite the system (fTSj) . defining the 

optimal vectors p, p', p", in the form: 

2(1 - pfWTx < l l Wl - p l Wp - 2(1 - p)*Wv, I G L, 

2(1 - p'YWTx < l l Wl - (p')Wp' - 2(1 - pO'W, Z G L', (29) 
2(Z - p'JWT* < l f Wl - (p"fWp" - 2(1 - p'JWv", I G L". 
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Note that we have suppressed the dependence of the sets L, L', L" on x , Xq, Xq. 
In fact, considering L, for example, it is clear that each candidate / G L must 
have sum of coordinates s and, in view of Lemma 14. 3\ its norm is bounded 
above by A + yK. Thus, taking 



K 



L = {I E Z K : / <A + VK, ^Z i = s }, 

i=i 

K K 

L' = {leZ K : \\l\\ <A + VK, J^a + l-^aJ}, 

i=i i=i 

K K 

L" = {I E Z K : / <A + VK, J2 l i = s + 2 ~12 a i}> 



i=l i=l 



(where I = (h,h, ■■■,^k))i we certainly do not miss any potentially optimal 
vectors p, p', p" by restricting the search to L, L', L", respectively. 

For each choice of rjx, rjx, s and of the vectors p, p', p", let P^Or ^e 
the set of all points in ^L_ VK satisfying (j^Hj) . Then 



K-l 3 3 

-Lff-1 



u j u u j u c p r. (3°) 

s=0 771=1 r\ K =l pGL p'GL' p"eL" 

where the sets on the right-hand side are disjoint (up to sets of a smaller 
dimension) . 

By (|T7|) and (J25j) . for all points in each subpolytope P^f v '^ we have the 
same optimal correction vectors 

yi = q-Tx, y 2 = q'-Tx, y 3 = q" - Tx, (31) 



where q = p — v, q' = p' — v', q" = p" — v". Hence, if the point x belongs to 
Pm^'r^x 1 then xn depends linearly on the coordinates {201}, {xq,k-i} ■ 

X n = l+(q - Tx) t H/(q - Tx) - 2(q' - Tx)*W{q - Tx) + 

+ (q ,/ -Tx) t H/(q"-Tx) (32) 
= 2(-q + 2q' - q")*WTx + canst, 
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where const = 1 + tfWq — 2(q')*H / q / + (q") t H / q".(Note that all coefficients 
on the right-hand side depend on 771, r/x, s, p, p', p".) 

We need to find the function F according to which the sequence (xn)n=i 
is asymptotically F-distributed. To simplify our notations, rewrite (|3U|) in 
the form 

i=i 

where each Pj is one of the polytopes P^^J 1 ' ■ Denote 

6 N = ({Nh}, {Nbjc-x}), N = 1, 2, .., 

and 

^ = {JV6N:^6Pi}, z = l,2,...,r. 

(The Aj's may intersect, as the Pj's may intersect on sets of a smaller di- 
mension. However, this will cause no problem as the intersections are sets 
of density in N. Alternatively, we may first decide in some arbitrary way 
where to place "problematic" integers.) 

Let 8$ be the subsequence of On, consisting of those elements On with 
N E Ai, 1 < i < r. That is, N ] is the iV-th element of (0n)n=i which 
belongs to A*. Let (Xn)n=x be the corresponding subsequence of (xn)n=x- 
If the point x lies in P, then it belongs to the subsequence (0^ )jv=i- Thus, 
by (J32j) there are exist affine functions ipi : R K_1 — > R such that 

X l 3 = &($})> l<«<r, iV = l,2,.... (33) 

Since the numbers bi,...,bx are linearly independent over Q, so are the 
numbers 1, b\, bx-i, and consequently the sequence (On) is uniformly dis- 
tributed modulo 1 in H K ~ 1 . Hence, each of the subsequences (0 N ^)^ =1 , 1 < 
i < r, is uniformly distributed in Pj. By (}3"3j) the sequence (xn)n=i * s the 
image of a uniformly distributed sequence in P; under the mapping ipi. Hence, 
letting (Yi, Y-i, Yk-x) be a (K — l)-dimensional random variable, uniformly 
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distributed in p, we see that (Xn )n=i is -^-distributed, where Pj is the dis- 
tribution function of 4>i(Yi, Y 2 , Yk-i)- 

According to jHl Thm 2.3], Pj is a piecewise polynomial function, each 
polynomial piece being of degree at most K — 1, and can be effectively com- 
puted. Since (0^)^=1 is uniformly distributed modulo 1 in R^ -1 , the density 
of each is the measure di of the set Pj. As Pj is a polytope, this measure 
can be effectively computed. By |2J Lemma 1], (xn)n=i is asymptotically 
P-distributed, where P = Y^i=i diF{. This completes the proof. 
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